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Abstract 

Multicriticality of the gonihedric model in 2 + 1 dimensions is investigated numerically. The 
gonihedric model is a fully frustrated Ising magnet with the finely tuned plaquette-type (four-body 
and plaquette-diagonal) interactions, which cancel out the domain-wall surface tension. Because 
the quantum-mechanical fluctuation along the imaginary-time direction is simply ferromagnetic, 
the criticality of the (2 + l)-dimensional gonihedric model should be an anisotropic one; that 
is, the respective critical indices of real-space (_L) and imaginary-time (||) sectors do not coin- 
cide. Extending the parameter space to control the domain-wall surface tension, we analyze the 
criticality in terms of the crossover (multicritical) scaling theory. By means of the numerical diag- 
onalization for the clusters with N < 28 spins, we obtained the correlation- length critical indices 
(^Xj^ll) = (0.45(10), 1.04(27)), and the crossover exponent <p = 0.7(2). Our results are comparable 
to {v±,v\\) = (0.482,1.230), and 4> = 0.688 obtained by Diehl and Shpot for the (d,m) = (3,2) 
Lifshitz point with the e-expansion method up to 0(e 2 ). 
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I. INTRODUCTION 



Recently, a thorough investigation of the Lifshitz point was made by Diehl and Shpot with 



2j; see also Refs. fl, fl, fl, |fi,E, ISj. The field theory 



the e-expansion method up to 0(e 2 ) 1 
for the Lifshitz point has an anisotropic dispersion like uj(k) = \ YlT=\^t + |X]f=m+i^ 2 > 
preventing us from going beyond order 0(e 2 ). Reflecting this anisotropy, the critical indices 
within the subspaces, i = 1,2, ... , m (_L) and i = m + 1, m + 2, . . . , d (||), are no longer 
identical. In Ref. the critical indices within each subspace are tabulated systematically 
for the generic values of (d,m). 

Such an anisotropic criticality is realized by the d- dimensional Ising model fully-frustrated 
within the m-dimensional subspace. The problem is that a naive computer simulation for 
the equilateral cluster does not yield adequate finite-size scaling. Rather, one has to adjust 
the shape of the cluster (that is, the system sizes of each subspace L\\ t ±) so as to fix the 
following scaled ratio to a constant value; 



const. 



(1) 



Here, the index z denotes the dynamical critical exponent, which characterizes the 
anisotropy. The significant point is that the exponent z itself is an unknown parameter, 
and it has to be determined through some preliminary analyses. After that, one is able to 
perform large-scale simulations. So far, the case of (d,m) = (3,1), namely, the axial-next- 
nearest-neighbor-Ising model, has been studied extensively by means of the Monte Carlo 



method 



0, Q 



111 ]. The simulation results are in agreement with the above-mentioned 
field-theoretical considerations as well as the series-expansion results Q, Q. 

In this paper, we consider the case of (d, m) = (3, 2). For that purpose, we investigate the 
ground-state phase transition of the gonihedric model in 2 + 1 dimensions. The gonihedric 
model is a fully frustrated Ising magnet with the finely tuned plaquette-type (four-body 
and plaquette-diagonal) interactions, for which the domain-wall surface tension vanishes; 
so far, the classical version has been studied in detail Q, 16, 17, 18]. Making a con- 



trast to the frustrated magnetism within the real space (_L), the quantum fluctuation along 
the imaginary-time direction (||) is simply ferromagnetic, and the ground-state criticality 
should be an anisotropic one. In Fig. [1] (a), we present a schematic phase diagram of the 
(2 + l)-dimensional gonihedric model subjected to the transverse magnetic field T and the 
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frustration j; we explain the details in Sec. [TTJ The multicritical point at j = 1, where the 
magnetism is fully frustrated, is our main concern. 

In order to simulate the (2 + l)-dimensional gonihedric model, we utilize the numerical- 
diagonalization method. This a ppr oach may have the following advantages. First, we 
implemented Novotny's method 19J to represent the Hamiltonian-matrix elements; this 
method is readily applicable to the quantum- mechanical system as well [20]. Owing to this 
method, we are able to treat an arbitrary number of spins iV = 8, 12, . . . , 28 constitut- 
ing the d = 2 cluster; note that conventionally, the number of spins is restricted within 
N(= L 2 ) = 9, 16, 25, .... Such an arbitrariness allows us to make a systematic finite-size 
scaling analysis. Second, the diagonalization method is free from the slowing-down problem; 
this problem becomes severe for such a frustrated magnetism, deteriorating the efficiency of 
the Monte Carlo sampling. Last, the constraint Lj_/L\\ — ► [Eq. (p0|] is always satisfied, 
because the system size along the imaginary-time direction is infinite L\\ — > oo; note that the 
system size along the imaginary time corresponds to the inverse temperature L\\ — 1/T — > oo. 

In fairness, it has to be mentioned that our research owes its basic idea to the following 
pioneering studies. First, an equivalence between the (2 + l)-dimensional fully frustrated 



magnetism and the (d, m 



21 



221. Second, in Ref. 



(3, 2) Lifshitz point was argued field-theoretically in Refs. 
23|, the biaxial-next-nearest-neighbor Ising model in d = 3 was 
studied with the Monte Carlo method. It was reported that the Lifshitz (multicritical) point 
collapses at zero temperature. On the contrary, the gonihedric model has an extra tunable 
parameter k. Setting k > 2, we attain desirable multicriticality as depicted in Fig. [T] (a). 

The rest of this paper is organized as follows. In Sec. HIl we explain the (2+l)-dimensional 
gonihedric model. To elucidate the underlying physics, we make an overview of the classical 
gonihedric model in d — 3. In Sec. IHIt we present the simulation results. The simulation 
scheme is explained in the Appendix. In Sec. IIV} we present a summary and discussions. 



II. QUANTUM GONIHEDRIC MODEL IN 2+1 DIMENSIONS: A REALIZATION 
OF THE (d,m) = (3,2) LIFSHITZ POINT 

In this section, we propose the (2 + l)-dimensional gonihedric model as a realization of 
the (d, m) = (3, 2) Lifshitz point. To elucidate the underlying physics, we make an overview 
of the original (classical) gonihedric model in d = 3. 
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A. Quantum gonihedric model in d = 2 



As mentioned in the Introduction, we propose the (2 + l)-dimensional gonihedric model as 
a realization of the (d, m) = (3, 2) Lifshitz point. To be specific, we consider the Hamiltonian 

h = -j 1 J2 °t°] - ^ E °M - J 3 E °t°K° z i - r E °f . ( 2 ) 

with the coupling constants Ji = k, J2 = —k/2, and J3 = (1 — k)/2. Here, the operators 
{erf} denote the Pauli matrices placed at the square-lattice points i. The summations XI (y) 5 
Xw/ij))) an d Sryjfc/i run over an possible nearest-neighbor, next-nearest-neighbor (plaquette 
diagonal), and plaquette-four-body spins, respectively. The transverse magnetic field T con- 
trols the amount of quantum fluctuations. At a certain T c , a ground-state phase transition 
may occur. 

As mentioned above, the gonihedric model has finely tuned coupling constants {Ji}, 
which cancel out the domain-wall surface tension. Actually, the domain-wall energy of the 
gonihedric model (apart from the off-diagonal term — T ^ erf) admits a geometric represen- 
tation E = ri2 + Anrii fl4| . Here, ri2 denotes the number of points where two domain walls 
meet at a right angle (domain- wall undulation), and is the number of points where four 
domain walls meet at a right angle (self-intersection point). That is, the parameter k con- 
trols the self-avoidance of the domain walls with the bending elasticity unchanged. (Notably 
enough, the interfacial energy lacks the surface-tension term. Accordingly, the domain-wall 
undulations are promoted, giving rise to a peculiar type of criticality.) The gonihedric model 
has a tunable parameter k with the zero surface tension maintained. This redundancy is 
an advantage over other frustrated magnetisms such as the biaxial-next-nearest-neighbor 
Ising model. We survey the regime k > 2, where we observed a clear indication of the 
Lifshitz-type criticality. 

In this paper, we extend the above-mentioned parameterization space. That is, introduc- 
ing a new controllable parameter j, we investigate the parameter space 

J\ = «, -h = — -, J3 = — ^— • (3) 

Note that at j = 1, the parameter space, Eq. ([2]), reduces to the above-mentioned one 
(original gonihedric model). Owing to the extension, the magnetic domain wall now acquires 
a finite domain- wall surface tension oc 1 — j. In other worlds, in terms of this extended 
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parameter space, we identify the Lifshitz point as a multicritical point; see the phase diagram 
in Fig. [I] (a). This viewpoint was proposed in Ref. [l6], where the authors investigate the 
criticality of the classical d = 3 gonihedric model with the cluster-variation method. In the 
next section, we will overview the properties of the classical gonihedric model, which may 
be relevant to the present study. 



B. Phase diagram of the classical gonihedric model: A brief overview 



Let us make an overview of the past studies of the (classical) gonihedric model. The model 
was introduced by Savvidy and Wegner as a lattice-regularized version of the string field 
theory [lj], [ijj] . However, recent developments dwell on the d = 3 case, aiming at a potential 
applicability to microemulsions. The criticality should belong to the Lifshitz point with 
the index (d, m) = (3,3), because the classical gonihedric model is isotropically frustrated. 
The (d, m) = (3, 3) criticality may be realized by the ternary mixture 24j] of water, oil and 



surfactant 



25 



26 



27 



28] ; actually, a crossover from the d = 3-Ising universality to an exotic 



one was reported in Refs. 



29 



30]. 



We present a schematic phase diagram of the (classical) d = 3 gonihedric model in Fig 
[U (b) la, |3l|. The Hamiltonian of the classical d = 3 gonihedric model is given by 

1 - K 



H 



(ij) (m 



An 



SiSjSkSi. 



(4) 



[ijki] 



(The Ising-spin variables {Si} are placed at the d = 3 lattice points.) We notice that 
the phase diagram resembles that of the quantum gonihedric model; the discrepancy as to 
j «-> —j is merely due to the difference of parameterization, and the subspace —j = 1/4 
corresponds to the fully-frustrated gonihedric model. 

A few remarks on the phase diagram follow: First, the Lifshitz point at — j = 1/4 is 



identified as an end-point of the critical branch (— j < 1/4) belonging to the d = 3 



sing 



universality. In fact, the multicritical (crossover) scaling theory applies successfully [16|, [31] 
to clarifying the nature of the Lifshitz point. (Direct numerical simulation at — j = 1/4 
appears to be rather problematic |32j.) We will accept this cross-over viewpoint as for 
the quantum gonihedric model. Second, in Refs. [17], [l8|, it was reported that for small 



k < 0.5, the multicritical point becomes a discontinuous one, accompanied with pronounced 



hysteresis. In particular, at K — 0, the model reduces to the so-called p-spin model 



33|, 



which is notorious for its slow relaxation to the thermal equilibrium (metastability). We 
found that a similar difficulty arises in the quantum gonihedric model. Hence, we devote 
ourselves to the large- k regime such as k > 2, where we observed a clear indication of the 
Lifshitz-type criticality. Last, the phase boundary separating the lamellar and ferromagnetic 
phases is (almost) vertical. This feature ensures that the multicritical point is located at 
—j = 1/4. The quantum gonihedric model possesses this property as shown in the next 
section. Actually, this is the most significant benefit of the parameterization scheme, Eq. 



III. NUMERICAL RESULTS 



In this section, we present the numerical results. Our aim is to estimate the critical 
indices (u±,uu) and <fi. As mentioned in the Introduction, we utilize Novotny's method 
to diagonalize the Hamiltonian (T5]) numerically. We explain the technical details in the 
Appendix. By means of this method, we simulated finite clusters with N < 28 spins. The 
linear dimension of the cluster L is given by the formula 

L = VN, (5) 

because the N spins constitute a d — 2 cluster. 



A. Finite-size scaling of the critical branch: d = 3-Ising universality 

In this section, we survey the critical branch j < 1; see Fig. [1] (a). We show that the 
criticality belongs to the ordinary d = 3-Ising universality class. This finding provides a 
foundation for the subsequent analyses with the crossover-scaling theory. 

In Fig. [21 we plot the Roomany-Wyld approximate beta function |34j 



rw (t) = 1 + HAE N (T)/AE N ^(T))/ ln(v^V/y^4) 
N ^d r AE N (T)d r AE N ^(T)/AE N (T)/AE N ^(Ty 

with the excitation energy gap AEn(T) for the system size N. Here, we fixed the self- 
avoidance parameter k — 2, and varied the frustration as j = —1.2, —0.8, —0.4, 0, and 
0.4. The zero point of the beta function indicates the location of the critical point r c (j). 
Basically, the critical branch depicted in Fig. [1] (a) follows from this analysis; afterward, we 
determine the critical point T c more precisely. 
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The slope of the beta function at T = T c yields an estimate for the inverse of the 
correlation-length critical exponent, 1/v. In Fig. [2], as a reference, we presented a slope 



(dotted line) 1/v = 1.5868 35] corresponding to the d = 3-Ising universality class. We see 
that the criticality is maintained in the d = 3-Ising universality class for a wide range of 
j. Actually, we obtained 1/v = 1.55, 1.54, 1.51, 1.46, and 1.39 for j = -1.2, -0.8, -0.4, 
0, and 0.4, respectively. These results demonstrate that the critical branch belongs to the 
d = 3-Ising universality class. It would be noteworthy that the shape of the beta function 
becomes distorted as j — > 1. That is, the regime exhibiting the slope 1/v shrinks gradually 
as j — > 1, indicating that a new type of criticality emerges at the multicritical point j = 1. 
Actually, we consider this crossover behavior rather in detail in the following sections. 
In Fig. [31 we present the approximate critical point T C (L 1 ,L 2 ) for [2/(Li + L 2 )] 1 ^ +UJ 



3J| with k = 2, j = -0.4, and 8 < Ni < N 2 < 28 (L lj2 = ^/N^)] here, we used the 
corrections-to-scaling exponent u = 0.821 and the exponent v = 1.5868~ 1 reported in Ref. 



35[. The approximate critical point T c (Li,L2) is determined by the zero point of the beta 



function. That is, it satisfies the equation 

V^A£ta(r c (Li,L2)) = ^ 2 AE N2 (T C (L U L 2 )). (7) 

From the least-squares fit to the data in Fig. [31 we obtained the critical point T c = 7.073(55) 
in the thermodynamic limit L — > oo. We make use of r c in the following scaling analyses. 



B. End-point singularity of the critical amplitude 

The above analysis indicates that the multicriticality at j = 1 is merely an end-point 
singularity of the ordinary d = 3-Ising critical branch. That is, the crossover-scaling theory 
should apply to clarifying the nature of the multicritical point. 

In this section, we consider the singularity of the critical amplitude of AE beside the 
multicritical point. The amplitude is defined by the relation 

AE^G ± U)\r-r c (j)\\ (8) 

The amplitude exhibits the singularity 

G±(A) ~ A^ir 1 ^, (9) 



7 



with the crossover exponent <fi. (As mentioned in the Introduction, the exponent u\\ denotes 
the critical index along the imaginary-time direction.) The variable A stands for the distance 
from the multicritical point 

A = l-j. (10) 

Here, we postulated that the multicritical point locates at j = 1, and we justify this claim 
in Sec. IIIID1 The above formula is a straightforward consequence of the crossover-scaling 
hypothesis 

ae« |r-r c p/(A/|r-r c |*). (11) 

Actually, this relation provides a definition of the crossover exponent <fi. 

To begin with, we determine the critical amplitude G + . In Fig. HI we plot the scaled 
energy gap (r - T^L^-AE /\T - T c \ u for k = 2, j = -0.4, and N = 8, 12, ... , 28. The 
critical point T c = 7.073 is determined in the above section, and likewise, we postulated the 



d = 3-Ising universality class \ jv = 1.5868 35J. The data collapse into a scaling-function 



curve. We again confirm that the phase transition belongs to the d = 3-Ising universality 
class. From the limiting value of the high-r side of the scaling function, we estimate the 
critical amplitude as G + = 4.28(8); here, we read off the value around the scaling regime 
(T — Y^L 1 ^ = 15, and accepted the data scatter among N = 20, 24 and 28 as an error 
indicator. 

Similarly, we determined G + for various values of j and k — 2, 4, and 6. In Fig. [5l we 
plotted the amplitude G + for A(= 1 — j) with the logarithmic scale. [In the cases of k — 2, 
4, and 6, we read off G + from the scaling plot at the scaling regime (r - T c )L l / v = 15, 
40, and 60, respectively. In the case of k — 2, we omitted the data of iV = 16 for its 
rather insystematic behavior particularly for small A.] In the plot, we also presented a slope 
(dotted line) of G + oc A 6 . We observe a signature of the power-law singularity with the 
exponent {y\\ — v)j<\> « 0.6 as A — > 0. Hence, we confirm that the crossover behavior (Q is 
realized in the vicinity of the multicritical point. In fact, from v = 0.63020(12) [35] and the 
present results, Eqs. (TI9]) and (TIBI) . obtained in Sec. IIII CI we arrive at the slope 

^-^= 0.59(42), (12) 



fairly consistent with the above observation. With use of G + calculated in this section, we 
crosscheck the validity of the critical indices obtained in the following section. 
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C. Finite-size-scaling analysis of (vx,vu) and <j) 

In this section, we make an analysis of each critical exponent with use of the crossover 
scaling, Eq. fTTT]) . 

First, we consider the Binder parameter 

with the magnetization M = YliLi °t- [Note that the simulation was not done right at the 
Lifshitz point; we calculated the data in the vicinity of the Lifshitz point (crossover scaling). 
Hence, the ferromagnetic order parameter M is still of use in the data analysis.] The symbol 
(. . . ) denotes the expectation value at the ground state. According to the crossover- scaling 
theory, the Binder parameter obeys the formula 

U = U((r-r c )L 1/vx ,AL^ v ±). (14) 

(Here, we made use of the fact that the Binder parameter is scale-invariant at the critical 
point.) As noted in the Introduction, the index with the subscript _L denotes the critical 
exponent within the real space. In Fig. [6j we present the crossover-scaling plot, (T—T C )L 1 ^ U± - 
U, with k = 2 and AL^^ U± = 8. Here, we set the scaling parameters v±_ = 0.45 and <fi = 0.7, 
where we found the best data collapse. Surveying k = 4 and 6 as well, we arrive at the 
estimates 

u ± = 0.45(10), (15) 

and 

= 0.7(2). (16) 

Second, we consider the energy gap AE. The energy gap obeys the crossover- scaling 
relation 

AE = L- z g((T - T c )L 1/u ^, AL* /vx ), (17) 

with the dynamical critical exponent z. In Fig. [71 we present the crossover-scaling plot, 
(r - T c )L l / v ^-L z AE, with k = 2 and AL^ = 8. Here, we set z = 2.3, and the other 
scaling parameters are the same as those of Fig. [61 Surveying k = 4 and 6 as well, we 
estimate the critical index as 

z = 2.3(3). (18) 
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Through z = v\\/u±, the above results lead to 



u\\ = 1.04(27). 



(19) 



Let us address a remark. As mentioned in the above section, the indices, Eqs. (jT6l) and 
( I19p . are consistent with the end-point singularity of G + , indicating the self-consistency of 
the present analyses. 

D. Phase transition between the ferromagnetic and lamellar phases 

The above analysis stems from the proposition that the multicritical point locates at 
j = 1; in other words, the phase boundary separating the ferromagnetic and lamellar phases 
is (almost) vertical. In this section, we justify this proposition. (Actually, this feature was 
confirmed in the case of the classical d — 3 gonihedric model [161] .) 

In Fig. [BJ we plot the ground-state energy per unit cell E /N with the system sizes 
N = 8,12, ... ,28 for k = 2 and T = 0.6; namely, we surveyed the regime slightly below the 
multicritical point. We observe a distinct signature of the first-order phase transition around 
j ~ 1, where the slope of E /N changes rather abruptly (level crossing). The transition 
point seems to converge into the regime 0.9 < j c < 1 as iV — > oo. Noticeably enough, the 
transition point is close to j = 1. 

We argue this behavior more in detail: First, the data E /N in j < j c (ferromagnetic 
phase) appear to reach the thermodynamic limit, whereas in j > j c (lamellar phase), the 
plots are still scattered insystematically. Possibly, the incommensurability of the lamellar- 
type structure (periodicity of the domain walls) causes such an irregularity. Surveying the 
cases of k = 2, 4, and 6, we found that the data of iV = 8, 16, and 24 are rather robust 
against this incommensurability effect. Hence, we conclude that the transition point locates 
within 0.9 < j c < 1. Last, we found that such a slight deviation of j c from j = 1 is 
negligible in the sense that the influence is less than the error margins. In other worlds, the 
parameterization, Eq. (J3j), is sensible to explore the multicriticality in terms of the cross-over 
scaling; this point was noted in the case of the classical gonihedric model [if]]. 
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IV. SUMMARY AND DISCUSSIONS 



We investigated the criticality of the (2 + l)-dimensional gonihedric model, Eq. (j2J), with 
the extended parameter space, Eq. (j3J). This extended parameter space allows us to survey 
the criticality in terms of the crossover-scaling theory; see Fig. [TJ(a). We employed Novotny's 
method to diagonalize the Hamiltonian. With use of this method, we treated an arbitrary 
(integral) number of spins N = 8, 12, . . . , 28. Because the quantum-mechanical fluctuation 
along the imaginary time direction is ferromagnetic, the criticality of the quantum gonihedric 
model should be an anisotropic one accompanied with the dynamical critical exponent z(= 
7^ 1- Our estimates for the critical indices are (u±, = (0.45(10), 1.04(27)) and = 
0.7(2). We also confirmed that the estimates are consistent with the end-point singularity 
of the critical amplitude G + ; see Fig. [5] 

As mentioned in the Introduction, Diehl and Shpot made an analysis of the Lifshitz point 
with the e-expansion method up to 0(e 2 ). Their conclusions for (d,m) = (3,2) are 



They also provided the convergence-accelerated results with the [1/1] Pade method; 



Our simulation data support their claim. 

Lastly, let us make a few comments on the advantages of the diagonalization approach. 
First, the numerical diagonalization is free from the slowing-problem problem, which deteri- 
orates the efficiency of the Monte Carlo sampling for the frustrated magnetism. Second, we 
do not have to worry about the constraint (pQ). The constraint is always satisfied, because 
the system size along the imaginary-time direction is infinite. However, the diagonalization 
method suffers from the severe limitation as to the available system sizes. In this paper, 
we surmount this difficulty with the aide of Novotny's method, which allows us to treat a 
variety of system sizes N = 8, 12, . . . , 28 sufficient to manage systematic finite-size scaling. 
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(i/ ± , z/||) = (0.387, 0.795) and (p = 0.686. 



(20) 




(21) 
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APPENDIX A: CONSTRUCTION OF THE HAMILTONIAN-MATRIX ELE- 
MENTS: QUANTUM NOVOTNY'S METHOD 



In this Appendix, we explain the simulation scheme. As mentioned in the Introduction, 



we applied the Novotny method 19J to diagonalizing the Hamiltonian (j2j). Novotny's method 
allows us to construct the Hamiltonian-matrix elements systematically for the cluster with 
an arbitrary (integral) number of spins N = 8, 12, . . . , 28; note that conventionally, the 
number of spins is restricted within N(= L 2 ) = 9, 16, 25, .... Originally, Novotny's method 



was formulated for the classical Ising model (transfer- matrix formalism) 19]. In Ref. 20], 



it was extended to adopt the quantum- mechanical interaction (Hamiltonian formalism). 



Here, we follow the notation of Ref. 20], and make a slight extension to incorporate the 
plaquette-type interactions; see Eq. (1A5I) . 

Before we commence a detailed discussion, we explain the basic idea of Novotny's method. 
In Fig. [9l we present a schematic drawing of a finite-size cluster for the d = 2 gonihedric 
model, Eq. (j2J). As seen in the figure, the spins {<Ji} (i = 1, 2, . . . , N) constitute a d = 1- 
dimensional (zig-zag) structure. This feature is essential for us to construct the cluster with 
an arbitrary (integral) number of spins N. The dimensionality is lifted to d = 2 by the long- 
range interactions over the v^Vth-neighbor distances; owing to the long-range interaction, 
the N spins constitute a \/N x \/N rectangular network effectively. (The significant point 
is that the number y/~N is not necessarily an integral nor rational number.) 

Let us formulate the above idea explicitly. To begin with, we set up the Hilbert-space 
bases {|ci, o-i, ... , &n)} {p% = ±1) for the quantum spins {erf} (i — 1,2, ... , N). The bases 
diagonalize the operator of; namely, the relation 

^IW> = ^-IW>, (Ai) 

holds. 

We decompose the Hamiltonian into two components 

H = H {D \{X i })+H {0 \T). (A2) 

The component 7i^ D \{ Jj}) describes the exchange interactions, depending on the coupling 
constants {Ji}- On the other hand, the component Ti^°\T) originates from the single-spin 
term, which depends on the transverse magnetic field T. The former component is a diagonal 
matrix, whereas the latter is off-diagonal. 
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First, we consider the diagonal component 7i^ D \ We propose the following formula 

n {D) = ^(h(Vn) + h(-Vn)). (A3) 

Here, the component H(v) is a diagonal matrix, which describes the vth- neighbor interaction 
among the iV-spin alignment. The diagonal elements are given by 

= (WI^(«)IM> = ({ai}\TP°\W)- (A4) 
Here, the matrix T denotes the plaquette-type interaction between the arrays {<jj} and {tj}; 

N ( h \ 

k=i ^ ' 

(A5) 

The operator P denotes the translational operator, which satisfies P\{<Ji}) = here, 
we imposed the periodic-boundary condition. Note that the operator insertion of P v in Eq. 
31) introduces the long-range interaction over the vth-neighbor pairs. The denominator 2 



in Eqs. (1A3R and (1A5[) compensates the duplicated sum. 

Lastly, we consider the off-diagonal component TiS°'. The matrix element is given by 

K!un } = ({^W { ° ] \{n}). (A6) 

The expression is quite standard, because the component TC^ simply concerns the individual 
spins, and has nothing to do with the connectivity among them. 

The above formulas complete our basis to simulate the Hamiltonian (T5]) numerically. The 
results are shown in Sec. IIIII 



[1] H. W. Diehl and M. Shpot, Phys. Rev. B 62, 12338 (2000). 

[2] M. Shpot and H. W. Diehl, Nucl. Phys. B 612, 340 (2001). 

[3] H. W. Diehl and M. Shpot, Phys. Rev. B 68, 066401 (2003). 

[4] M. M. Leite, Phys. Rev. B 67, 104415 (2003). 

[5] M. M. Leite, Phys. Rev. B 68, 066402 (2003). 

[6] R. M. Hornreich, M. Luban, and S. Shtrikman, Phys. Rev. Lett. 35, 1678 (1975). 

[7] C. Mergulhao, Jr and C. E. I. Carneiro, Phys. Rev. B 58, 6047 (1998). 

[8] C. Mergulhao, Jr and C. E. I. Carneiro, Phys. Rev. B 59, 13954 (1999). 

13 



[9] W. Selke, Z. Phys. B 29, 133 (1978). 
[10] K. Kaski and W. Selke, Phys. Rev. B 31, 3128 (1985). 
[11] M. Pleimling and M. Henkel, Phys. Rev. Lett. 87, 125702 (2001). 
[12] J. Oitmaa, J. Phys. A: Math. Gen. 18, 365 (1985). 
[13] Z. Mo and M. Ferer, Phys. Rev. B 43, 10890 (1991). 
[14] G.K. Savvidy and F.J. Wegner, Nucl. Phys. B 413, 605 (1994). 

[15] R.V. Ambartzumian, G.S. Sukiasian, G.K. Savvidy, and K.G. Savvidy, Phys. Lett. B 275, 99 
(1992). 

[16] E.N.M. Cirillo, G. Gonnella, and A Pelizzola, Phys. Rev. E 55, R17 (1997). 

[17] G. Koutsoumbas and G.K. Savvidy, Mod. Phys. Lett. A 17, (2002) 751. 

[18] D. Espriu, M. Baig, D.A. Johnston, and Ranasinghe P.K.C. Malmini, J. Phys. A: Math. Gen. 

30, 405 (1997). 
[19] M.A. Novotny, J. Appl. Phys. 67, 5448 (1990). 
[20] Y. Nishiyama, Phys. Rev. E 75, 011106 (2007). 

[21] A. Dutta, B. K. Chakrabarti, and J. K. Bhattacharjee, Phys. Rev. B 55, 5619 (1997). 
[22] A. Dutta, J. K. Bhattacharjee, and B. K. Chakrabarti, Eur. Phys. J. B 3, 97 (1998). 
[23] B. Pal and S. Dasgupta, Z. Phys. B 78, 489 (1990). 

[24] K.A. Dawson, M.D. Lipkin, and B. Widom, J. Chem. Phys. 88, 5149 (1988). 
[25] J. S. Huang and M. W. Kim, Phys. Rev. Lett. 47, 1462 (1981). 
[26] P. Honorat, D. Roux, and A. M. Bellocq, J. Phys. (Paris) Lett. 45, L-961 (1984). 
[27] R. Aschauer and D. Beysens, Phys. Rev. E 47, 1850 (1993). 

[28] H. Seto, D Schwahn, M Nagao, E Yokoi, S Komura, M Imai, and K Mortensen, Phys. Rev. E 
54, 629 (1996). 

[29] D. Schwahn, K. Mortensen, H. Frielinghaus, and K. Almdal, Phys. Rev. Lett. 82, 5056 (1999). 
[30] P. Stepanek, T. L. Morkved, K. Krishnan, T. P. Lodge, and F. S. Bates, Physica A 314, 411 
(2002). 

[31] Y. Nishiyama, Phys. Rev. E 70, 026120 (2004). 

[32] R.K. Hellmann, A.M. Ferrenberg, D.P. Landau, and R.W. Gerling, Computer Simulation Stud- 
ies in Condensed- Matter Physics IV, edited by D.P. Landau, K.K. Mon, and H.-B. Schiittler 
(Springer- Verlag, Berlin, 1993). 

[33] J.D. Shore and J.P. Sethna, Phys. Rev. B 43, 3782 (1991). 

14 



r t 

K=2 (T=0) 




FIG. 1: (a) A schematic drawing of the ground-state phase diagram of the (2 + l)-dimensional 
gonihedric model, Eq. ([2]), with the self- avoidance parameter k = 2. We aim to investigate the 
multicritical point at j = 1. (b) As a comparison, we present the phase diagram of the d = 3 
(classical) gonihedric model, Eq. with k = 1 Q^lsil; here, the parameter T denotes the 

temperature. The phase diagram is essentially the same as that of the quantum-mechanical model; 
the discrepancy j <-> — j is due to the difference of parameterization. 

[34] H.H. Roomany and H.W. Wyld, Phys. Rev. D 21, 3341 (1980). 
[35] Y. Deng and H.W.J. Blote, Phys. Rev. E 68, 036125 (2003). 
[36] K. Binder, Z. Phys. B: Condens. Matter 43, 119 (1981). 



15 



2 4 6 8 10 12 

r 



FIG. 2: The beta function (3^ (T) © is plotted for k = 2 and various j. The symbols, +, x, *, 
□ , and ■, denote the data for j = 0.4, 0, —0.4, —0.8, and —1.2, respectively. For a comparison, 
we presented a slope (dotted line) corresponding to the d = 3-Ising universality class (y = 0.6294 



35]). We see that the criticality is maintained to be that of the d = 3-Ising universality class for 



a wide range of j; see text for details. 
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FIG. 3: The approximate critical point T C (L 1 ,L 2 ) is plotted for [2/(Li + L 2 )] l/u+u with 8 < N x < 
N2 < 28, k = 2, and j = —0.4; the corrections-to-scaling exponent u) = 0.821 and the exponent 
1/u = 1.5868 are taken from Ref. [35J. The least-squares fit to these data yields T c = 7.073(55) in 



the thermodynamic limit. 
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FIG. 4: The scaling plot for the energy gap, (T — T^L 1 ^- AE/\T — T c \ u , is shown. The parameters 



are the same as those of Fig. [3l We postulated the d = 3-Ising universality class v = 0.6302 
The symbols, +, x, *, □, ■. and o, denote the system sizes of N = 8, 12, 16, 20, 24, and 
28, respectively. We confirm that the transition belongs to the d = 3-Ising universality class. 
Furthermore, from the plateau in the high-r side, we obtain an estimate for the critical amplitude 
G + = 4.28(8); see text for details. 
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FIG. 5: The correlation-length critical amplitude G + ([9j) is plotted for various A(= 1 — j) and 
k = 2, 4, and 6. The symbols x, □, and o denote the data for k = 2, 4, and 6, respectively. 
As a reference, we presented a slope (dotted line) of G + oc A 0,6 . The data indicate a power-law 
singularity, Eq. Q, with the exponent — v)/<j)Xi 0.6; see text for details. 
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FIG. 6: The crossover-scaling plot (O, (r - T c )L l / v ±-U, for k = 2 and AL^ = 8 is shown. 
Here we set v± = 0.45 and <j> = 0.7, for which we found the best data collapse. The symbols, +, 
x, *, □, ■. and o, denote the system sizes of N = 8, 12, 16, 20, 24, and 28, respectively. 
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FIG. 7: The crossover-scaling plot (fTT|) . (r - Y c )L l / VjL -L z AE, is shown. Here, we set z = 2.3, and 
the other scaling parameters are the same as those of Fig. [6j The symbols, +, x, *, □, ■, and o, 
denote the system sizes of N = 8, 12, 16, 20, 24, and 28, respectively. 
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FIG. 8: The ground-state energy per unit cell Eq/N is plotted for k = 2 and V = 0.6. The symbols, 
+, x, *, □, ■. and o, denote the system sizes of N = 8, 12, 16, 20, 24, and 28, respectively. There 
occurs a transition separating the ferromagnetic and lamellar phases. The transition point seems 
to converge into the regime 0.9 < j c < 1 as N — > oo; see text for details, 
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FIG. 9: Construction of the spin cluster for the gonihedric model ([2]). As indicated above, the 
spins constitute a d = 1-dimensional alignment (i = 1,2,... , N), and the dimensionality is 

lifted to d = 2 by introducing the bridges (long-range interactions) over the (iV 1 / 2 )th -neighbor 
pairs. 
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